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Abstract. Parameter dependent non-Hermit ian quantum systems typically not 
only possess eigenvalue degeneracies, but also degeneracies of the corresponding 
eigenfunctions at exceptional points. While the effect of two coalescing 
eigenfunctions on cyclic parameter variation is well investigated, little attention 
has hitherto been paid to the effect of more than two coalescing eigenfunctions. 
Here a characterisation of behaviours of symmetric Hamiltonians with three 
coalescing eigenfunctions is presented, using perturbation theory for non- 
Hermitian operators. Two main types of parameter perturbations need to be 
distinguished, which lead to characteristic eigenvalue and eigenvector patterns 
under cyclic variation. A physical system is introduced for which both behaviours 
might be experimentally accessible. 
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1. Introduction 



The coalescence of eigenvalues and the corresponding eigenvectors of operators 
at so-called exceptional points (EPs) in parameter space plays a crucial role in 
various branches of physics, reaching from classical mechanics [T]|2] and optics [^11 



12 



13 . Since Hermit ian operators have a complete set 



to quantum mechanics 
of eigenfunctions, EPs are a phenomenon that can only occur for non-Hermitian 
operators. In quantum mechanics, however, the Hamiltonian is traditionally 
demanded to be Hermitian for the description of closed systems. Non-Hermiticity, 
and with it EPs, can enter quantum mechanics in different ways: Most importantly, 
as an effective description of open systems 12 , or in the form of PT-symmetric 
quantum mechanics 14 aiming at a generalised description of closed systems flSj. 
Further, the analogy between the Schrodinger and the Helmholtz equations allows 
beautiful experimental realisations of non-Hermitian quantum systems in microwave 



cavities 16 -19 , and the fact that the spatial propagation of light in wave guide arrays 



is described by an equation similar to the time-dependent Schrodinger equation led 
to experimental realisations of non-Hermitian quantum dynamics in optical systems 
20)[2l]. 

EPs have been investigated in one form or another for decades |3l^-24 



and 

they have attracted considerable attention in recent years, with the increasing interest 
in non-Hermitian quantum theories. The behaviour of eigenvalues and eigenvectors 
under cyclic parameter variation in the presence of EPs has been investigated in 
pO], a nd the predictions have been confirmed experimentally 
16!T7lfT9l 



detail theoretically 125 



in microwave cavities 



The influence of EPs on dynamics [3l pT|[3T||33] , 

leading to typical quadratic behaviour for two coalescing eigenfunctions could also 
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be demonstrated in experiments pS). The experimental studies on EPs were so far 
confined to systems which are governed by wave equations that ahow to mimic the 
quantum behaviour, such as the above mentioned microwave cavities, laser cavities, or 
optical wave guide structures. Only recently there have been suggestions for possible 



experimental observations of EPs in genuine quantum systems 30 33 - 35 



While EPs at which two eigenfunctions coalesce are well investigated, little 
attention has hitherto been paid to the coalescence of more than two eigenfunctions. 
Among the few exceptions are 3T|[36 39 . The typical effect on the dynamics at an EP 
of more coalescing eigenfunctions is a straight forward generalisation of the behaviour 
at a two fold degeneracy [31^, the influence on the cyclic behaviour, on the contrary. 



is more elaborate. In 38 an interesting chiral behaviour of the eigenfunctions in the 
neighbourhood of three coalescing eigenfunctions has been reported. The topological 
structure of three interacting eigenfunctions in the form of neighbouring EPs of each 



two eigenfunctions has been investigated in 40 . In this context it was found that the 



interchange of the eigenfunctions when encircling several EPs can be quite different 
from the expected behaviour around a degeneracy of three eigenfunctions, which is 
often identified with a third root branch point in the literature 1370381. However, while 



the third root behaviour is typical^ in the sense specified below, it is by no means the 
only possible structure [22|[39 |[4l|[42 . Thus, although the coalescence of more than two 
eigenfunctions is less commonly found in physical systems than the standard EP due to 
the higher codimension ^ , the nontrivial features associated with such degeneracies 
make it desireable to gain a deeper understanding into the related structures. The 
purpose of the present paper is to initiate these investigations. 

Here we present a study of the quantum signatures in the neighbourhood of three 
coalescing eigenvalues for symmetric Hamiltonians. We focus in particular on the 
behaviour of eigenvalues and eigenvectors under parameter variation, and the resulting 
interchanges between eigenfunctions and possible additional geometric phases. Using 
generalised perturbation theory for non-Hermitian operators we show that there are 
two distinct fundamental scenarios. We suggest an experimental setup, consisting of 
three coupled wave guides, in which both behaviours could in principle be observable. 

The paper is organised as follows: We begin with a short overview on the 
mathematical definition of exceptional points and the related Jordan block theory 
in Section |2] In Section |3] we review the general series expansion of eigenvalues 
and eigenvectors in dependence on a complex perturbation parameter. We analyse 
the resulting patterns of eigenvalues under variation of parameters, in particular, 
for the cyclic variation of a complex parameter. We then turn to the behaviour of 
the eigenvectors under cyclic parameter variation, and the corresponding geometric 
phases in Section [4] We present a physical example in Section [5] In this context 
we draw attention to possible problems when relying on adiabatic assumptions in 
experiments [43||45] . In Section|6]we briefiy comment on a more general two-parameter 
perturbation. We summarise our results in Section [7| 



2. Exceptional Points and Jordan block theory 

Let us begin with a brief summary of the mathematics of exceptional points, where 
we also introduce the notations used here. For details on Jordan chains and related 
topics see, e.g., 1 . For convenience we work in a finite-dimensional Hilbert space. 
Hamiltonians acting on this space are represented by complex non-Hermitian matrices 
depending on a number of system parameters. An exceptional point (EP) is a 
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point in parameter space at which at least two eigenvalues and the corresponding 
eigenvectors of the Hamiltonian coalesce. That is, at these points the Hamiltonian 
is not diagonalisable, but it can be brought to Jordan normal form via a similarity 
transformation H = RJR~^^ where J is block diagonal, consisting of Jordan blocks 
Jn associated with the eigenvalues A^: 



1 



A, 



1 

An / 



(1) 



If N eigenvalues and their corresponding eigenvectors coalesce while the remaining 
eigenvalues are nondegenerate, we refer to the configuration as an EPA^. At an EPA^ 
the Jordan normal form of H contains one Jordan block of length N and other than 
that only blocks of size one, corresponding to the nondegenerate eigenvectors. 

If an eigenvalue is A/'-fold degenerate with only one linearly independent 
eigenvector, there exist — 1 associated eigenvectors that form the so-called Jordan 
chain: 

(F-A)|mo)=0, {H-\)\uj) = \uj_^), (2) 

where {H — \)\x) = \un-i) is insoluble. 

When dealing with non-Hermitian Hamiltonians one further has to distinguish 
left and right eigenvectors, where the left eigenvectors are commonly defined as the 
right eigenvectors of the adjoint operator corresponding to the complex conjugate 
eigenvalue. For a nondegenerate spectrum the left and right eigenvectors form a 
biorthogonal basis. In the degenerate case in addition to the eigenvectors the left and 
right Jordan chains are necessary to provide a basis for the state space. In analogy 
with the left eigenvector the left Jordan chain is defined as the Jordan chain of the 
adjoint operator W: 

{H^ -X)\vo) = 0, (Ft-A)|«,.) = |«,_i), (3) 

where {W — X)\x) = I'^at-i) is insoluble. 

In the present study we confine ourselves to the investigation of symmetric 
matrices as a first step towards the understanding of the physics of three coalescing 
eigenvectors. Note, however, that the generalisation to nonsymmetric matrices can 
lead to qualitatively new effects, in particular, with respect to the geometric phase, 
as has been investigated in detail for EP2s 19,29 . In the case of symmetric matrices 
the left eigenvectors and Jordan chain vectors reduce to the complex conjugates of the 
right eigenvectors and Jordan chain: \vj) = \uj). Thus the scalar product between left 
and right eigenvectors reduces to the standard inner product of real Hilbert spaces, 

{Vj\Uk 



ujuk^ which in the present context is often referred to as the c-product 



12 



For convenience we will from now on omit the bracket notation and denote the right 
eigenvectors and Jordan chains simply by uj. At an EP it holds that uo"^uo = 0, 



which can be interpreted as self- orthogonality 12 



The Jordan chain vectors are a priory not uniquely defined, however, with 
additional normalisation conditions the ambiguity can be removed (up to a sign) [l]. 
In addition to the self-orthogonality condition they fulfil a number of further relations, 
depending on their length. In what follows we shall mainly be concerned with EP2s 
and EP3s, thus we will confine ourselves to the discussion of Jordan chains of length 
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two and three here. For a Jordan chain of length two we have the self-orthogonahty 
condition 

EP2 : uo^uo = 0. (4) 

Further, it is convenient to impose the conditions 

EP2 : uo^ui = 1 and ui^ui = 0, (5) 

which uniquely determine the Jordan chain (up to a sign). For an EPS, that is, a 
Jordan chain of length three, these relations change and one finds instead 

EP3 : uo^uo = 0, uq^ui = 0, and uo^U2 = ui^ui. (6) 

To remove the ambiguity the following relations are imposed: 

EP3 : uo^U2 = 1, and U2^ui = = U2^U2. (7) 

As we shall see below the structure of the Jordan chain at the EP is the key for the 
topological phases related to encircling an EP. The typical interchange of eigenvalues 
and eigenfunctions, which is often discussed as part of the same phenomenon, is in 
fact due to the perturbative behaviour around the EP, which we shall discuss in the 
following section. 



3. Perturbation of parameters around exceptional points 

One of the most characteristic features of EPs is the behaviour of the eigenvalues 
under a perturbation in the neighbourhood of the EP. For Hermitian operators the 
perturbed eigenvalues can be expanded into a Taylor series in the perturbation 
parameter. This behaviour holds also for non-Hermitian operators, as long as 
the perturbation is performed around a nondegenerate eigenvalue or a degenerate 
eigenvalue corresponding to nondegenerate eigenvectors [22^. If, however, a system 
is perturbed around an EP the series expansion is in general a Puiseux series that 
might involve fractional exponents. Note, however, that contrary to common folklore, 
this does not necessarily have to be the case. Instead, the following theorem on the 
"typical behaviour" of eigenvalues under a perturbation around an EP holds ^ : 

Theorem 1. Let A be a linear perturbation of Ao: 
A=Ao^zAl, (8) 

where z G C, and Aq is a full Jordan block of size N with the eigenvalue Aq. // the 
perturbation matrix Ai = (akj) (the matrix representation of Ai in the Jordan basis 
of Aq ) is such that a^i ^ 0^ then the set of perturbed eigenvalues of A consists of the 
N branches of an analytic function, i.e. 

oo 

A,(z)=Ao + ^A^,z^/^, (9) 

where Xia- = cil/i^ , and where a labels the different branches of the N-th root. 

The condition a^i 7^ is crucial here. Although it is often argued that "most 
perturbations fulfil this condition" , it turns out that this is not true for many physical 
situations, as for example in lattice models described by band matrices. If, however. 
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Figure 1. Typical cases for the two possible perturbation scenarios of the 
eigenvalues around an EPS. Shown here are the real (left column) and the 
imaginary (right column) parts of the eigenvalues of a 3 x 3 matix as a function of 
a complex perturbation parameter in the neighbourhood of an EP3. The upper 
panel shows the familiar third-root behaviour. The lower panels show an example 
for which a square-root perturbation is observed for two of the eigenvalues, 
whereas the third eigenvalue is constant. This behaviour is usually connected 
to EP2s. The examples depicted here correpond to the Hamiltonian ( |27| with 
-y = I = V and a = b = 2c (top), and a = —b = 26 (bottom). 



this condition is not fulfilled, the perturbed eigenvalues can split into different cycles 
each of which consists of the branches of an analytic function of the form 

oo 

A,(^) = Ao + ^A^,W^ (10) 



where k < and the different values of k sum up to 22 41 . In the case of 
symmetric matrices, a similar expansion is valid for the eigenvectors. 

One might argue that for an EP2 this is of little consequence, as the only 
two possibilities in this case are a square-root expansion and a conventional Taylor 
expansion, where the latter can be viewed as a special case of the former for which 
all odd coefficients vanish. For degeneracies of more than two eigenvalues and 
eigenfunctions, however, the situation is quite different. In particular, it follows 
that for the case of an EPS there are two possible power expansions (apart from 
the Taylor expansion, which can be viewed as a special case of one of the other 
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Figure 2. Third-root behaviour of eigenvalues in dependence on a perturbation 
parameter around an EP3, and resulting eigenvalue patterns under a cyclic 
parameter variation. The left and the middle columns show the real and the 
imaginary parts, respectively, as functions of the complex perturbation parameter. 
The solid black line illustrates the eigenvalue curves that are observed when the 
parameter is varied along a one-dimensional closed loop. The right column shows 
the corresponding eigenvalue trajectories in the complex plane. Three distinct 
behaviours are observed for a loop outside the EP (top), passing through the EP 
(middle), and encircling the EP (bottom). 

expansions) for the eigenvalues in the perturbation parameter z. The first (better 
known) expansion consists of a single cycle with the three perturbed eigenvalues being 
the three branches of a third root expansion. In the second (less appreciated) scenario 
one of the eigenvalues has a conventional Taylor expansion, and the remaining two 
are given by the branches of a square-root expansion. This second scenario is similar 
to the characteristic behaviour of an EP2. Typical examples of the two cases are 
depicted in figure [l] in the top and bottom panel, respectively. The figure shows the 
real and the imaginary parts of the eigenvalues of a 3 x 3 matrix perturbed around an 
EP3 as a function of the complex perturbation parameter. The situation in the upper 
panel corresponds to a third-root behaviour, while the situation in the lower panel 
corresponds to the second scenario. Note that in the latter case one of the eigenvalues 
is independent of the perturbation parameter up to the leading order of the other two 
eigenvalues. 

This perturbative behaviour gives rise to typical patterns of eigenvalues for cyclic 
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Figure 3. As in figure [2] however, for a square-root dependence on the 
perturbation parameter. 

parameter variation. It is well known that for an EP2, due to the square-root 
behaviour, the two eigenvalues interchange when the EP is surrounded and form a 
single closed loop in the complex plane, while they describe two seperate closed curves 
if the EP is not enclosed. In the case of an EPS the typical patterns of eigenvalues 
under parameter variation for the two possible scenarios are depicted in figures [2] and 
[Sj respectively. When encircling an EPS along a parameter loop that corresponds to 
a triple-root perturbation, the three eigenvalues interchange, which results in a single 
closed curve formed by the three eigenvalues in the complex plane depicted in the 
lower right panel of figure [2] This is the standard behaviour usually connected to an 
EPS. However, the second perturbation scenario yields eigenvalue patterns commonly 
connected with EP2s, as depicted for an example in figure [Sj 

The explicit perturbative expressions for the eigenvalues and eigenvectors can be 
obtained from a generalised perturbation method 1 , by inserting an ansatz for the 
series expansion into the eigenvalue equation and equating powers of the perturbation 
paramter z. This yields the following leading order behaviour for an EP2: 

A = Ao + Aiz^/'+o(z^/'), with Ai = (uo^HiUo)'/', (11) 
u = uo + AiUi^i/2+ 0(^1/2). (12) 

That is, the leading order correction for the eigenvectors is collinear to the first Jordan 
vector ui at the EP. This is in fact similar for an EPS, where in the case of a third 
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root series expansion one finds: 

A = Ao + Ai^i/^ + o(zi/^), with Ai = (uo^HlUo)'/^ (13) 
u = Uo + Aiuiz^/^ + o(z^/^). (14) 

For both cases the perturbation series has to be modified in the case uo"^^iUo = 0, 
which in fact corresponds to the condition a^i in Theorem [l] For an EP2 the 
alternative behvaiour is given by a Taylor expansion in the perturbation parameter [l] , 
which we shall not discuss here in further detail. 

In the case of an EPS the alternative behaviour is the square-root expansion of 
two of the eigenvalues and eigenvectors, and a Taylor expansion for the remaining 
eigenvalue and eigenvector, which explicitly yields: 

' Ao ± A^'^zi/2 + 0(^1/2)^ with A^'^ = (2uiTHiUo)^/2 



A 



Ao + Af^z + o(zi/^), with X\ 



l/2\ 



,(2) 



UQ^HiG-^HiUQ 
2uoTHiUi 



uo + A^'^uizV2 
uo + (Af^ui 



o{z) 



■o{z), 



(15) 



(16) 



where the invert ible operator G = i^o — '^ol ~uoU2"^ has been introduced by removing 
the kernel of the noninvertible operator — AqI. 

The fact that the leading order perturbations of the eigenvectors are branches of 
an analytic function leads to an interchange of the eigenvectors under cyclic parameter 
variation. While the square-root expansion of the second EPS scenario leads to an 
interchange between two eigenvalues and the two corresponding eigenvectors after one 
cycle just as in the case of an EP2, the phase that is acquired by the eigenvectors after 
two cycles, is different in these two cases, due to the different Jordan block structure 
at the EP. In what follows we shall discuss this behaviour in more detail. 



4. Behaviour of the eigenvectors under cychc parameter variations 

In Hermitian quantum mechanics, a system initially in an eigenstate returns to this 
initial state after an adiabatic cyclic variation of parameters around an eigenvalue 
degeneracy. During this cyclic variation it acquires an additional geometric phase; the 
celebrated Berry-phase 46 . In the presence of EPs, however, the situation is more 
intricate. It is a well known result, which we shall briefiy review in what follows, that 
surrounding an EP2 in parameter space, the two eigenvectors interchange after one 
cycle and two cycles are needed to return to the initial eigenstate. In addition the 



states pick up a phase of tt after two cycles 16j|25 29 . Another important difference 
to the Hermitian case has been pointed out recently in 4S -45 : In the presence of non- 
Hermiticities the adiabatic theorem might break down in the sense that non-adiabatic 
corrections do not have to vanish. As a consequence fundamental difficulties for the 
observation of the phases or even the interchange behaviour using adiabatic parameter 
variation are to be expected. Here we shall for the moment not be concerned with 
this issue, and first analyse the situation for abstract cyclic parameter variations. 

In this section we shall investigate the effect of a cyclic variation of parameters 
around an EPS. In particular, we consider the variation of one complex parameter 
around a loop in the complex plane encircling the EPS. We find that the geometric 
phase is determined only by the order of the EP, while the number of parameter cycles 



Signatures of three coalescing eigenf unctions 



9 



that are necessary for the eigenvectors to return to their initial state is governed by the 
Puiseux expansion form of the eigenvectors, and can thus depend on the perturbation 
type around an EPS. The situation for more general parameter perturbations is briefly 
discussed in Section |6l 

To begin with, let us review the situation for EP2s. The leading order of the 
eigenvectors in the neighborhood of the EP in dependence on the complex parameter 
z, which we shall decompose into its amplitude and phase, z = re^"^, has the form 

u± =uo± V^ae^^/^ui, (17) 

where a G C is constant. 

To determine the behaviour under a cyclic variation of the parameter z it is 
necessary to fix the normalisation of the eigenvectors consistently for all parameters. 
In particular to obtain meaningful results for the phase, care has to be taken that the 
parallel transport condition (that the derivative of the eigenfunction with respect to 
the parameter be orthogonal to the eigenfunction itself) is fulfilled. Here we use the 
standard normalisation for the eigenvectors of non-Hermitian operators: u^u± = 1, 
which also automatically ensures the parallel transport condition. Using the self- 
orthogonality condition Q and the Jordan chain relations ([5| we find in leading order 
of r that 

u±^u± = ±2av^e^^/^ (18) 
and thus we get a suitable expression for u±: 

u+ = . (19) 

When the parameter performs one cycle, that is, when the angle traverses the interval 
[0,27r] , the numerators in (19) interchange. This phenomenon is purely due to the 
square root perturbation behaviour of the eigenvectors. From this it follows that up 
to a phase, which is determined by the denominators, the eigenvectors return to their 
initial state after two cycles of parameters in the case of an EP2. We will see shortly 
that a similar behaviour can prevail for an EPS. Now we turn to the behaviour of the 
denominator in ( 19 ). Recalling that the complex square root displays a half line branch 
cut in the complex plane, one finds that either the positive or the negative argument 
in the square root passes through the branch cut. This induces a minus sign to one 
of the denominators when going through one cycle, and a minus sign to the other 
when passing through the second cycle. This corresponds to the geometric phase of tt 
that the eigenvectors acquire after two cycles. In summary, combining the behaviours 
of both the denominator and the numerator we obtain the well-known result for the 
behaviour of the eigenvectors under successive cyclic parameter variations (denoted 
by O): 

u+ O u_ O -u+ O -u_ O u+. (20) 

The two eigenvectors interchange for each cycle, where one picks up an additional 
minus sign. Thus, after two cycles the system returns to its initial eigenvector up to a 
geometric phase of tt. This phase is intrinsically related to the Jordan chain structure 
at the exceptional point. 

We now proceed in an analogous manner to determine the behaviour in the 
presence of EPSs. We begin with the standard triple root case (14), where the three 
relevant eigenvectors in leading order in the parameter z = re^^ have the form 

u^'') = uo + ae"=2V3^i/3 gW3„^_ (21) 
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The superscript (k) denotes the different branches of the complex third root. Choosing 
again the standard normaUsation (u^^^)-^u^^^ = 1, and using the Jordan chain 
properties ([6| and ^ , we find the fohowing expression for the normahsed eigenvectors 
in leading order of r: 

We observe that the numerators are cyclically interchanged in each cycle. Thus, after 
three cycles, the eigenvectors return to themselves, up to a possible phase, which is 



determined by the behaviour of the denominator in (22). The crucial part in the 



denominator is the factor Ve^^^^/^ g2i</)/3^ When cf) describes a full cycle, each of the 
arguments in the square root describes two thirds of a full circle in the complex plane. 
As a consequence, in each cycle two of the arguments pass through the branch cut, 
and thus, the denominators cyclically interchange and two of the three pick up a 
minus sign. If we choose the signs of the eigenvectors in the initial configuration 
appropriately, this minus sign can be absorbed in the definition of the eigenvectors, 
and we find the following pattern for successive parameter cycles: 

u(i) 0u(2) Ou(3) Ou(i). (23) 

The three eigenvectors cyclically interchange for each parameter cycle and thus return 
to their initial state after three cycles, where none of them picks up a geometric phase. 
While the three cycles required to recover the original configuration are related to the 
triple root structure of the parameter dependence, the fact that no phase is acquired 
is due to the Jordan block structure at the EPS, and we will find the same behaviour 



for the second perturbation scenario (16). In this case we first consider the two 



eigenvectors with the square root behaviour, which take the same form as in the EP2 
case: 

u± = uo iav^e^^/^ui. (24) 



However, due to the different Jordan chain relations at the EPS, instead of (19) we 
find for the normalised eigenvectors in the lowest order of r: 

V a r e^^ 

As in the case of the EP2, the numerators interchange in each parameter cycle. The 
behaviour of the normalisation, however, differs from the EP2 case. For each cycle of 



paramaters the argument in the square root of the denominator in (25) performs a 
whole cycle in the complex plane, and hence the denominator picks up a minus sign. 
Again we can absorb the sign by an appropriate choice of the initial eigenvectors. 
Putting these together, we find for the behaviour of two of the eigenvectors under 
consecutive parameter cycles: 

u+ O u_ O u+ (26) 

The two eigenvectors interchange for each cycle of parameters and thus return to their 
initial state after two cycles just as in the EP2 case, since this interchange is encoded in 
the Puiseux expansion. However, in contrast to the EP2 case, none of the eigenvectors 
picks up a geometric phase. A study of the remaining eigenvector in a similar way 
shows, that it returns to itself after each cycle while picking up a minus sign, that is, 
a phase of tt. 
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In summary, we found that the form of the Puiseux expansion of the eigenvectors 
determines the number of parameter cycles that are necessary for the eigenvectors 
to return to their initial state, whereas the order N of the degeneracy determines 
the geometric phase. Note that the results here are based on the perturbation 
around the degeneracy of three eigenfunctions, which is not restricted by the actual 
dimensionality of the considered Hilbert space. However, modifications for larger 
circles in parameter space, for which higher orders of the perturbation expansion 



become relevant are possible 28 . Finally, the observed phenomena should directly 



generalise to degeneracies of more than three eigenvectors. 

5. A physical example: EP3s in three coupled wave guides 

Let us now present an experimental setup in the context of optics, for which the two 
different behaviours of an EP3 might be observable. In particular, we consider light 
propagation in three approximately parallel coupled wave guides. In the paraxial 
approximation the spatial light propagation is described by an equation analogous to 
the time dependent Schrodinger equation, where one spatial direction takes the role of 
the time, and the potential is replaced by the refractive index of the wave guides, which 
can be complex in the case of absorbing or amplifying structures. Recently an EP2 



has been experimentally observed in a similar structure of two wave guides 20,21 



For more details on the analogy between quantum dynamics and optical structures. 



see, e.g., |6l^ 49 



The three wave guide system we consider here is constructed to resemble a 
quantum system with a Hamiltonian of the form 

/ a-2i7 V^v \ 
H= V2v V2v , (27) 
\ V^v 6 + 2i7 / 

with 7,'u G R, and a,6 G C. This corresponds to a situation where each of the wave 
guides supports a single mode. The real and the imaginary parts of the diagonal 
elements are related to the refractive index and the absorption and amplification in 
the wave guides, and can thus be changed by choosing different materials, or varying 
the amplification strength or depths of the wave guides (within a range for which 
the single mode approximation remains valid). Neighbouring wave guides are coupled 
via the evanescent field between them. Thus, the coupling parameter v is related to 
the distance of the wave guides (the coupling is large if the wave guides are close 
to each other and approaches zero, if they are far apart). The numerical values 
appearing in the Hamiltonian have been chosen for convenience such that the system 
displays an EP3 for the values 7 = v and a = = b. Since the time dependence in 
the Schrodinger equation corresponds to a spatial variation in the optical context, 
an adiabatic variation of parameters can be achieved by a variation of the wave 
guide parameters along the beam propagation direction on a length scale that is long 
compared to the intrinsic scales of the system. The parameters necessary for an 
observation of the predicted phenomena should be within experimental reach, since 
the loop surrounding the exceptional point, and thus the corresponding parameter 
variations, can be arbitrarily small. Furthermore, the presence of the EPS relies only 
on the ratio of absorption and coupling, and the differences in the refractive indices 
of the three waveguides, but not on absolute values. In the context of Hermitian 
quantum dynamics many comparative models with adiabatic time dependence have 
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0.05 -005 



Figure 4. Real (left) and imaginary (right) parts of the eigenvalues of the 
Hamiltonian ( |27| for j = 1 = v in dependence on the parameters a and b. The 
solid line shows the eigenvalues along a curve in parameter space that encircles 
the exceptional point. The figure on the right shows the corresponding eigenvalue 
trajectory. 



been realised in wave guide structures [iT. However, it has recently been shown [4^ 45 
that the non-Hermiticity, that is, the presence of loss and gain in the wave guides, 
can destroy the expected adiabaticity. This might prevent a measurement of the 
interchange behaviour and the phases under adiabatic parameter variation. Instead 
one might have to be contended with the measurement of instantaneous eigenvalues 
and modes. 

The two different scenarios discussed in the present paper can in principle be 



observed in a system described by the Hamiltonian (27). It can easily be verified that 



the eigenvalues and eigenvectors of (27) around the EPS have the typical third root 
dependence on the complex parameter z if we choose a = b = whereas they show 
the square root behaviour for the choice a = —b = z. In fact the two different cases 
depicted in figure [l] are the eigenvalues of the present example, where j = v = 1. 



6. Outlook: More general two-parameter perturbations 

The results presented here rely on the Puiseux expansion in the neighbourhood of 
degeneracies. It is interesting to note that in contrast to Taylor expansions around 
nondegenerate eigenvalues, there are mathematical subtleties leading to qualitatively 
new behaviour for multi-paramter perturbations around Jordan blocks. This is 
connected to the fact that there is no direct generalisation of Puiseux expansions to 
several variables (see, e.g., 50 ). Thus, if parameters are varied in a way which cannot 
be combined into a single complex parameter, the behaviour cannot be captured by 
the perturbative considerations presented above. A detailed analysis of the possible 
structures for these cases, however, goes beyond the scope of the present paper, and is 
an interesting topic for future investigations. Here, we shall only present an example 
to illustrate some structures that might be expected. 



For this purpose we return to the model (27), and consider a variation of the 
independent variables a, 6 G R, which cannot be combined into a single effective 
complex parameter. In figure [4] we show the resulting eigenvalues for 7 = 1 = v. 
It can be seen that although the eigenvalues form a three-sheeted Riemann surface, 
we do not find the typical structures discussed above of either a triple root or a 
square root and an additional constant eigenvalue. Instead, the eigenvalues form 
a more intricate structure with multiple additional connections, which resembles a 
square-root behaviour in dependence on the parameter b in the plane a = and a 
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triple-root behaviour in dependence on the parameter a in the plane b = 0. The 
figure also shows the resulting eigenvalue pattern for a cyclic parameter variation 
(a, 6) = r(cos((/)), sin(0)), with a small constant r G H and G [0,27r]. Although 
the EP is surrounded, the resulting eigenvalue trajectories show a self- crossing. This 
is usually related to parameter loops crossing an EP, whereas here it results from 
the additional connections in the imaginary parts of the eigenvalues. The related 
behaviour of the eigenvectors is correspondingly more elaborate and will be discussed 
in a future study. 

7. Discussion 

In summary, we have presented the behaviour of eigenvalues and eigenvectors of 
symmetric non-Hermitian operators in dependence on a complex parameter in the 
neighbourhood of exceptional points at which three eigenvalues and the corresponding 
eigenfunctions coalesce. In particular, we have studied the patterns resulting from 
cyclic parameter variations. We have shown that there are two basic scenarios, related 
to the two possible Puiseux expansions of the eigenvalues and eigenvectors. They lead 
to characteristic interchanges of eigenvalues and eigenfunctions for consecutive cycles. 
In addition, there can be geometric phases, which are related to the degeneracy itself 
and independent of the particular parameter variation. We have presented a possible 
experimental setup in which both behaviours might be observable. Finally, we have 
illustrated that the patterns resulting from general cyclic variation of parameters 
(which cannot be expressed as the variation of one effective complex parameter) 
can show additional qualitatively new features. Further interesting generalisations 
are expected for non-symmetric operators. An extension of the present study to 
exceptional points of higher orders, along the lines presented here, is a promising 
topic for future investigations. 
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